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ABSTRACT 

We report the detection of a planetary companion with a minimum mass of 
msini = 0.0771 M Jup = 24.5 M e to the nearby (d = 9.4 pc) M2.5V star GJ 176. 
The star was observed as part of our M dwarf planet search at the Hobby- Eberly 
Telescope (HET). The detection is based on 5 years of high-precision differen- 
tial radial velocity (RV) measurements using the High-Resolution-Spectrograph 
(HRS). The orbital period of the planet is 10.24 d. GJ 176 thus joins the small 
(but increasing) sample of M dwarfs hosting short-periodic planets with mini- 
mum masses in the Neptune-mass range. Low mass planets could be relatively 
common around M dwarfs and the current detections might represent the tip of 
a rocky planet population. 

Subject headings: planetary system — stars: individual (GJ 176) — techniques: 
radial velocities 



1 Based on observations obtained with the Hobby-Eberly Telescope, which is a joint project of the Uni- 
versity of Texas at Austin, the Pennsylvania State University, Stanford University, Ludwig-Maximilians- 
Universitat Miinchen, and Georg-August-Universitat Gottingen. 
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1. Introduction 

Over the past few years our preliminary knowledge on the frequency of extrasolar planets 
in the low mass part of the HR-diagram has increased significantly. An ever increasing 
number of M dwarfs is being monitored by various Doppler searches with high radial velocity 
(RV) precision. This has led to several discoveries of M dwarf planets that cover an enormous 
and surprising range in mass, almost comparable to the mass range of the planets in our 
Solar System. Jovian planets were found around GJ 876 (Delfosse et al. 1998, Marcy et 
al. 1998, 2001), GJ 849 (Butler et al. 2006) and GJ 317 (Johnson et al. 2007), and Neptune- 
mass planets around GJ 436 (Butler et al. 2004), GJ 581 (Bonfils et al. 2005b) and GJ 674 
(Bonfils et al. 2007). The low primary masses of M dwarfs combined with state-of-the-art 
RV precision even allowed the detection of additional planets with minimum masses below 
10 M ffi (so-called "Super-Earths") in the GJ 876 (Rivera et al. 2005) and GJ 581 (Udry et 
al. 2007) systems. 

In this paper we report the detection of a new low mass planet in a 10.24 d orbit around 
the M dwarf GJ 176. This is the fourth planet in the class of Neptune-mass companions 
orbiting M dwarfs. GJ 176 is already the third (previously unknown) planet host in our 
HET M dwarf planet search (Endl et al. 2003, 2006) sample. 

2. Stellar parameters of GJ 176 

GJ 176 (HD 285968, HIP 21932, LHS 196) is a V = 9.97 M2.5Ve star at a distance of 
9.4 parsec, according to the Hipparcos parallax of 106 ±2.5 mas (Perryman et al. 1997). The 
star has a B — V of 1.51 and an absolute V magnitude of 10.08. The 2MASS magnitudes for 
GJ 176 are J = 6.462, H = 5.824 and K = 5.607 mag (Cutri et al. 2003). Using the V-band 
and K-band mass-luminosity relationships of Delfosse et al. (2000), we estimate a mass of 
0.48 and 0.50 M , respectively. We adopt the mean value of 0.49 ± 0.014 M as the mass of 
the star. 

A photometric study by Weis (1994) didn't find significant variability for GJ 176. The V- 
band scatter measured by Weis (1994) is 0.006 mags, equal to the measurement uncertainties. 
The ROSAT All-Sky-Survey catalog of nearby stars (Hiinsch et al. 1999) reports a moderate 
coronal X-ray emission level of 3 x 10 27 ergs _1 . GJ 176 is thus a moderately active star, 
possibly exhibiting starspots and flares, a quite typical behavior for M dwarfs. Rauscher 
& Marcy (2006) measured the equivalent widths of the Ca II H and K lines with 0.97 ± 
0.10 A and 0.69 ± 0.08 A, respectively. Based on the V-band mass-metallicity relationship 
of Bonfils et al. (2005a) we estimate an [Fe/H]= —0.1 ± 0.2 for GJ 176, so roughly solar 
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metallicity. 

3. Observations and RV results 

We observed GJ 176 as part of our on-going Doppler search for planets around M dwarfs 
(Endl et al. 2003, 2006) using the HET (Ramsey et al. 1998) and its HRS spectrograph 
(Tull 1998). We started observations of GJ 176 in 2003 October 15 and collected a total of 
28 spectra over 5 years. All observations were performed using our standard planet search 
setup and data reduction pipeline decribed in detail in Cochran et al. (2004). We use the 
common I 2 cell technique to obtain high precision differential RV measurements (e.g. Butler 
et al. 1995, Endl et al. 2000). 

Fig. [1] shows the time series of our HET/HRS RV measurements with the small secular 
acceleration of the RV of 0.36 ms^yr -1 , as computed from the Hipparcos parallax and 
proper motion information, already subtracted. The data have an overall rms-scatter of 
9.84 ms _1 and average internal errors of 4.69 ± 0.63 ms _1 . The total scatter is more than 
twice the measurement uncertainties, indicative of intrinsic RV variability of this target. The 
HET RV data are listed in Table [TJ 

4. Period search and orbital solution 

We searched the RV data of GJ 176 for any significant periodicities using the classic 
Lomb-Scargle periodogram (Lomb 1976 ; Scargle 1982). Fig. [2] displays the resulting power 
spectrum. A strong peak is visible at a period of 10.24 d. We estimate the false-alarm- 
probability (FAP) of this signal with the bootstrap randomization scheme (e.g. Kiirster et 
al. 1997). After 100,000 bootstrap re-shuffling runs we find that the FAP of the 10.24 d 
signal is only 0.0004. 

As the next step we use Gaussfit (Jeffereys et al. 1988) to find a Keplerian orbital solution 
to our RV data. A circular orbit fit yields a \ 2 of 35.2 (x 2 ed = 1.47) and a residual RV scatter 
of 5.57 ms _1 . A slightly better fit is obtained with an eccentric orbit with e = 0.23 ± 0.13: 
X 2 of 32.86 (x? e d = 1-49) and residual rms of 5.32 ms _1 . The large uncertainity in e and 
the fact that a circular orbit yields a lower x 2 ed urges us to remain cautious concerning the 
reality of the non-circular orbit. However, future observations will allow us to determine 
whether the orbit is indeed eccentric. A moderately eccentric orbit could be an indication 
for additional planets around this star and thus warrants intensive follow-up monitoring. 

Fig. [3] shows the RV measurements and both orbital solutions (circular and eccentric) 
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phased to the best-fit period. The circular orbital solution yields a minimum mass for the 
companion of 0.077 ± 0.012 Mj up = 24.5 ± 3.9 M®, while an eccentric orbit would lower the 
minimum mass slightly to 0.076 ± 0.010 Mj up = 24.1 ± 3.1 M®. The orbital parameters are 
summarized in Table [2J 

4.1. Stellar activity versus planet hypothesis 

The case of GJ 674 has demonstrated clearly that star spots can introduce low level 
signals in high precision RV data of M dwarfs (Bonfils et al. 2007). Can the GJ 176 signal 
also be caused by rotational modulation due to star spots and not a planet? This scenario is 
very unlikely because a rotation period of 10.2 d would mean that GJ 176 would rotate more 
than 3 times faster than GJ 674, and this should lead to a higher activity level than GJ 674. 
This is not the case, as GJ 176 has a slightly lower coronal X-ray emission than GJ 674 
(Hiinsch et al. 1999). As mention before, Weis (1994) didn't detect photometric variability 
in GJ 176 and we also did not find any significant variability or periodicity in the Hipparcos 
photometry for this star (the highest peak at 3.7 d has a FAP of 1.5%.) 

Moreover, the vast majority of the data of Bonfils et al. (2007) was collected over a 
relatively short period of time (~ 200 d), where rotational modulation can mimic a Keplerian 
signal of a planet because the active regions on the star remain constant for this period of 
time. However, over a larger amount of time, active regions will reconfigure and emerge at 
different stellar longitudes causing a phase as well as amplitude change in the RV signal. 
The GJ 176 signal remains stable in phase and amplitude over 5 years, which makes the 
planet hypothesis much more plausible as the origin of the RV signal. 

5. Discussion 

Remarkably, GJ 176 is already the third M dwarf with a Neptune-mass companion that 
was included in our HET sample of 60 M dwarfs. GJ 436 and GJ 581 were part of our 
M dwarf sample before the planets around them were announced. The frequency of short- 
periodic Neptune-mass planets around M dwarfs in our HET sample is hence 5%, which 
is higher than the frequency of hot Jupiters around FGK-type stars of ~ 1.2% (Marcy et 
al. 2005). Neptunes and Super-Earths could be relatively common around low mass stars. 
Of course, with M dwarfs we are still limited by small number statistics, as compared to the 
few thousands of FGK-type stars already observed by various Doppler surveys around the 
world. 
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We are currently increasing the size of the HET M dwarf sample, and combined with the 
results of other programs observing M dwarfs, we should be able to derive important con- 
straints for planet formation models. The present information about Neptune-mass planets 
on short-period orbits suggests that they may be the tip of the terrestrial planet distribution, 
for several reasons. First, a number of hot Neptunes have 2 or 3 sibling gas giant planets 
orbiting at much greater distances: p 1 Cnc, \i Ara, and GJ 876 all have such systems. This 
planetary system architecture is the same as our own solar system, with inner terrestrial 
planets and outer gas giant planets, and suggests that the hot Neptunes in these systems 
formed inside their gas giants, making them likely to be rocky planets. Theoretical models 
of the collisional accumulation of terrestrial planets predict that rocky planets as massive as 
about 21 M e would result from a protoplanetary disk with a surface density of gas and solids 
high enough to form a gas giant planet by core accretion in a few million years (Wether- 
ill 1996; Inaba et al. 2003). Formation of the Neptunes interior to their gas giants seems to 
be a much more likely scenario than formation as ice giants on orbits outside the orbits of the 
gas giants, followed by migration somehow past the gas giants to their current short-period 
orbits. While GJ 876 is the only M dwarf of these three stars, one might expect a common 
explanation for its hot Neptune as a rocky world, as well as those of GJ 436 and GJ 581. 
This interpretation bodes well for the eventual detection of Earth-mass planets on habitable 
orbits around low mass stars. For M dwarfs with masses of < 0.2 M Q we already have the 
sensitivity to detect msinz 1 M e planets in the habitable zone (Endl & Kurster 2007, in 
prep.). 

While M dwarfs appear to have hot Neptunes at a higher rate than G dwarfs, they appear 
to have short-period gas giant planets at a somewhat smaller rate than G dwarfs (Endl et 
al. 2006 ; Johnson et al. 2007). Similarly, microlensing surveys have detected planetary 
companions at asteroidal distances to M dwarf stars at frequencies that suggest that planets 
at such locations are more likely to be Neptune-mass than Jupiter-mass. These results 
may simply be a result of M dwarfs have lower mass protoplanetary disks than G dwarfs, 
with the consequent result that the planets that form tend to have a mass distribution 
that is shifted downward. However, this apparent preference for cold Super-Earths over 
Jupiters can also be explained at present by both the core accretion and disk instability 
formation mechanism for gas giant planets. In the case of core accretion, the lengthened 
orbital periods around M dwarfs could lead to most cores growing too slowly to accrete 
significant gaseous envelopes, making failed cores the usual outcome, rather than gas giants 
(Laughlin et al. 2004). Disk instability explains the appearance of both gas giants and Super- 
Earths at asteroidal distances around M dwarfs by appealing to rapid formation of gaseous 
protoplanets by disk instability, followed by conversion to ice giants (cold Super-Earths) by 
photoevaporation at asteroidal distances in region of high mass star formation (Boss 2006). 
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Because most M dwarfs form in regions of high mass formation, they should be orbited at 
asteroidal distance primarily by Super-Earths. M dwarfs that form in region of low mass star 
will have gas giant planets at those distances, if disk instability is operative (Boss 2006). It 
remains for these theoretical speculations to be tested by completing the extrasolar planetary 
census. 

This material is based on work supported by the National Aeronautics and Space Admin- 
istration under Grants NNG04G141G, NNG05G107G issued through the Terrestrial Planet 
Finder Foundation Science program and Grant NNX07AL70G issued through the Origins of 
Solar Systems Program. The Hobby-Eberly Telescope (HET) is a joint project of the Uni- 
versity of Texas at Austin, the Pennsylvania State University, Stanford University, Ludwig- 
Maximilians-Universitat Miinchen, and Georg-August-Universitat Gottingen. The HET is 
named in honor of its principal benefactors, William P. Hobby and Robert E. Eberly. We 
would like to thank the McDonald Observatory TAC for generous allocation of observing 
time. 
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Fig. 1. — 5 years of HET/HRS RV measurements of GJ 176. The data have a scatter of 
9.84 ms" 1 and an average measurement uncertainty of 4.69 ms _1 . 
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Fig. 2. — Lomb-Scargle periodogram of the RV data of GJ 176. The highest peak is at 
10.24 days and has a false-alarm-probability of 0.0004. The lower panel displays the window 
function of our observations. 
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Table I. Differential radial velocities for GJ 176 from the HET/HRS 



JD [2,400,000+] dRV [ms^ 1 ] t [m s ^ ] .ID [2,400,000+] dHVIms" 1 ] t [m s ^ ] 
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Fig. 3. — The best-fit orbital solutions are shown as solid line (circular) and dashed line 
(eccentric) along with the HET/HRS RV measurements phased to the orbital period of 
10.24 d (the data are plotted twice for a second cycle). The semi-amplitude K of 11.6 ms _1 
corresponds to a minimum mass of 24.5 M ffi for the companion. All parameters of the 
orbital solution are summarized in Table [2j The residual scatter around the fit is 5.57 ms -1 
(circular) and 5.32 ms" 1 (eccentric orbit). 
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Table 2. Parameters of GJ 176 b for circular (left column) and eccentric orbit (right 

column) 



Parameter Value Value 



P 10.2369±0.0039 days 10.2366±0.0038 days 
T 2454550.6672±0.39 BJD 2455037. 7979±1.1 BJD 

K ll.62dzl.61 ins" 1 ll.72il.62ms- 1 

e 0.0±0.0 (fixed) 0.232±0.127 

lo O.OiO.O degrees (fixed) 210.4±32.5 degrees 

Msinz 0.0771±0.0122 M Jup 0.0757±0.0096 M Jup 

Msini 24.5±3.9 M e 24.1±3.1 M© 

a 0.0727±0.0007 AU 0.0727±0.0007 AU 

RMS 5.57 ms- 1 5.32 ms- 1 



